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The generation of a momentum-dependent gluon mass proceeds through a sophisticated im- 
plementation, at the level of the Schwinger-Dyson equation for the gluon propagator, of the 
Schwinger mechanism, whose central dynamical ingredient is the nonperturbative formation of 
longitudinally coupled massless bound-state excitations. In addition to triggering the aforemen- 
tioned mechanism, these excitations introduce poles in the various off-shell Green's functions of 
the theory, in such a way as to maintain the Slavnov-Taylor identities intact in the presence of mas- 
sive gluon propagators, acting effectively as composite Nambu-Goldstone bosons. In this work 
we focus on the dynamics leading to the actual formation of such bound states. Specifically, we 
derive and solve numerically an approximate version of the homogeneous Bethe-Salpeter equa- 
tion governing the wave function of this special bound state. It is found that this integral equation 
admits physically meaningful non-trivial solutions, indicating that the QCD dynamics produce 
one of the crucial ingredients required for the gauge-invariant generation of a gluon mass. 
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1. Introduction 



It is by now a well-established fact that large-volume lattice simulations in the Landau gauge 
yield a gluon propagator that reaches a finite non-vanishing value in the deep infrared [|l], ^, ^, 
0, HI |6|, 0]. Without a doubt, the most physical way of explaining this observed finiteness is to 
invoke the mechanism of dynamical gluon mass generation, first introduced in the seminal work 
of Cornwall ||8|], and subsequently studied in a series of articles |^, 10, |ll|]. In this picture the 
fundamental Lagrangian of the Yang-Mills theory (or that of QCD) remains unaltered, and the 
generation of the gluon mass takes place dynamically, through the well-known Schwinger mecha- 
nism 1 12, 13, 14, |T5[ |l^, [l^, 18], without violating any of the underlying symmetries (for related 
contributions and alternative approaches, see, e.g., [19, 2C, 21 , 2^ 23, 24, 25, 26]). 



The way how the Schwinger mechanism generates a mass for the gauge boson (gluon) can 
be seen most directly at the level of its inverse propagator, A^^(^^) = q^[\ + /n(^^)], where n(^) 
is the dimensionless vacuum polarization. According to Schwinger's fundamental observation, 
if n(^^) develops a pole at zero momentum transfer {q^ = 0), then the vector meson acquires a 
mass, even if the gauge symmetry forbids a mass term at the level of the fundamental Lagrangian. 
Indeed, if n(^^) = /q^, then (in Euclidean space) A~^(^^) = q^ + m^, and so the vector meson 
becomes massive, A^'(O) = m^, even though it is massless in the absence of interactions (g = 0, 

n = o) plpl]. 



The key assumption when invoking the Schwinger mechanism in Yang-Mills theories, such 
as QCD, is that the required poles may be produced due to purely dynamical reasons; specifi- 
cally, one assumes that, for sufficiently strong binding, the mass of the appropriate bound state 
may be reduced to zero [0, [ll, |16[ l^, l^]. In addition to tri ggering the Schwinger mechanism, 
these massless composite excitations are crucial for preserving gauge invariance. Specifically, the 
presence of massless poles in the off-shell interaction vertices guarantees that the Ward identities 
(Wis) and Slavnov Taylor identities (STIs) of the theory maintain exactly the same form before 
and after mass generation (i.e. when the the massless propagators appearing in them are replaced 
by massive ones) [^, 17, 18, 11]. Thus, these excitations act like dynamical Nambu-Goldstone 



scalars, displaying, in fact, all their typical characteristics, such as masslessness, compositeness, 
and longitudinal coupling; note, however, that they differ from Nambu-Goldstone bosons as far as 
their origin is concerned, since they are not associated with the spontaneous breaking of any global 
symmetry Finally, every such Goldstone-like scalar, "absorbed" by a gluon in order to acquire 
a mass, is expected to actually cancel out of the S'-matrix against other massless poles or due to 
current conservation [14,[T5|, |l^, |17, 18]. 

The main purpose of this presentation is to report on recent work [^, where the central 
assumption of the dynamical scenario outlined above, namely the possibility of actual formation 
of such massless excitations, has been examined. Specifically, the entire mechanism of gluon 
mass generation hinges on the appearance of massless poles inside the nonperturbative three-gluon 
vertex, which enters in the Schwinger Dyson equation (SDE) governing the gluon propagator. 
These poles correspond to the propagator of the scalar massless excitation, and interact with a pair 
of gluons through a very characteristic proper vertex, which, of course, must be non vanishing, or 
else the entire construction is invalidated. The way to establish the existence of this latter vertex is 
by finding non-trivial solutions to the homogeneous Bethe-Salpeter equation (BSE) that it satisfies. 
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(ai) (as) 



Figure 1: The "one-loop dressed" gluon contribution to the PT-BFM gluon self-energy. White (black) 
circles denote fully dressed propagators (vertices); a gray circle attached to the external legs indicates that 
they are background gluons. Within the PT-BFM framework these two diagrams constitute a transverse 
subset of the full gluon SDE. 



2. Basic concepts 

The full gluon propagator A"^y{q) = S^^A^viq) in the Landau gauge is defined as 



where 



PMq)=snv — — 



(2.1) 



(2.2) 



is the usual transverse projector, and the scalar cofactor A{q^) is related to the (all-order) gluon 
self-energy n^v(^) = ^^^v(^)^(^^) through 

A-\q^)=q^ + in{q^). (2.3) 

One may define the dimensionless vacuum polarization 11 (^^) by setting Tl{q^) = q^Ti{q^) so that 

A-\q^)=q^[\+in{q^)]. 



(2.3) becomes 



Alternatively, one may define the gluon dressing function J{q^) as 

A-\q^)=q^J{q^). 



(2.4) 



(2.5) 



In the presence of a dynamically generated mass, the natural form of A ^{q^) is given by (Euclidean 
space) 

A-\q^)=q^J{q^)+m\q^), (2.6) 

where the first term corresponds to the "kinetic term", or "wave function" contribution, whereas 
the second is the (positive-definite) momentum-dependent mass. If one insists on maintaining the 
form of (2^) by explicitly factoring out a q^, then 



A-\q^)=q' 



Jiq') + 



m 



q^ 



(2.7) 



and the presence of the pole, with residue given by m^(0), becomes manifest. 
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The Schwinger mechanism is integrated into the SDE of the gluon propagator through the 
form of the three-gluon vertex. In particular, a crucial condition for the realization of the gluon 
mass generation scenario is the existence of a special vertex, to be denoted by Vap.v{(l,r,p), which 
must be completely longitudinally coupled, i.e. must satisfy 

P«'«(^)P^'^(r)P^'^(p)y„^v(g,r,p) =0. (2.8) 

The role of the vertex Va^v{^,r,p) is instrumental for maintaining gauge invariance, given 
that the massless poles that it must contain in order to trigger the Schwinger mechanism, act, at the 
same time, as composite, longitudinally coupled Nambu-Goldstone bosons. Specifically, in order 
to preserve the gauge invaiiance of the theory in the presence of masses, the vertex Va^v{Q,r,p) 
must be added to the conventional (fully-dressed) three-gluon vertex Jra^v{(l,r,p), giving rise to 
the new full vertex, Jr'fy,^y{q,r,p), defined as 

^'anv{^:r,p) = Taiiv{q,r,p)+Vaiiv{q,r,p) . (2.9) 

Gauge invariance remains intact because IT' satisfies the same WI (or STI) as IT before, but now 
replacing the gluon propagators appearing on their rhs by massive ones; schematically, A^^ — > A 



'm ' 



where the former denotes the propagator given in (2.5), while the latter that of (2.6). 



To see this in detail, let us employ the formalism provided by the synthesis of the pinch tech- 



nique (PT) |28|, |29p with the background field method (BFM) [|30[]. In this framework, the natural 
quantity to consider is the vertex BQQ, to be denoted by ^aiiviq^r^p), connecting a background 
gluon (B) with two quantum gluons {Q). With the Schwinger mechanism turned off, this vertex 
satisfies the WI 

^«ir«^v(^,r,/p) =/7(/)P;,v(/')-rV(r2)P^v(r), (2.10) 

when contracted with respect to the momentum of the background gluon. Then, gauge invariance 
requires that 

q"Va^,y{q,r,p) = m\r^)P^y{r) - m^{p^)P^y{p) , (2.11) 

SO that, after turning the Schwinger mechanism on, the corresponding WI satisfied by IT' would 
read 

q"^'a^viq^f^p) = q"[^{q,f,p)+v{q,r,p)]^^^ 

= [p'j{p^)-m\p')]P,y{p) - [r^J{r')-m\r')]P,y{r) 

= A;,\p^)P^y{p)-A^\r^)P^yir), (2.12) 

which is indeed the identity in Eq. ( [2. lO] ), with the aforementioned replacement A^^ A,^^ en- 
forced. The remaining STIs, triggered when contracting ir'^^^(^, r,p) with respect to the other two 
legs are realized in exactly the same fashion. 

The next step is to insert ir'„^^,(^, r,/?) into the SDE equation satisfied by the gluon propa- 
gator, see Fig. |l|. Then, a rather elaborate analysis [[TI|] gives rise to an integral equation for the 
momentum-dependent gluon mass, of the type 

m^{q^) = l^m^{k^)K{q,k), (2.13) 

where the kernel K survives the ^ ^ limit, i.e., lim^^o K{q,k) ^ 0, precisely because it includes 
the term l/q^ contained inside Va^v{q,r,p)- 
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3. Structure of the pole vertex 



rQ„„(9,r,p) = 





(«i) 



(«2) 



(«3) 



(04) 



Figure 2: The SDE for the BQQ vertex which connects a background gluon {B) with two quantum gluons 
(0. 

The main characteristic of the vertex V , which sharply differentiates it from ordinary vertex 
contributions, is that it contains massless poles, originating from the contributions of bound-state 
excitations. Specifically, all terms of the vertex V are proportional to 1/^^, 1/ r^, 1 / p^, and products 
thereof. Such dynamically generated poles are to be clearly distinguished from poles related to 
ordinary massless propagators, associated with elementary fields in the original Lagrangian. 

To see how such poles enter into the vertex, let us focus on the general structure of the SDE 
for the BQQ vertex. With the Schwinger mechanism turned off, the various multiparticle kernels 
appearing in this SDE have a complicated skeleton expansion (not shown here), but their common 
characteristic is that they are one-particle-irreducible with respect to cuts in the direction of the 
momentum q 

When the Schwinger mechanism is turned on, the structure of the kernels is modified by the 
presence of composite massless excitation, described by a propagator of the type i/q^, as shown in 
Fig. I The sum of such dynamical terms, coming from all multiparticle kernels, shown in Fig. ^, 



constitutes a characteristic part of the vertex V , to be denoted by JJ in Eq. ( |3.2[ ), namely the part 
that contains at least a massless propagator i/q^. The remaining parts, to be denoted by R, contain 
massless excitations in the other two channels, namely r^/r^ and py/p^ (but no qa/q^), and are 
not relevant for the purposes of this presentation. Thus, 



with 



Vanv{q,r,p) = Ua^iv{q,r,p)+Ranv{q,r,p), 
Uanviq, r,p) = qa {vigiiv + V2<?^j<?v + V^PiiPv + V4r^^v + ^sr^p^ 



(3.1) 



(3.2) 



where the Vi are form factors depending on the various momenta. 

At this point we can make the nonperturbative pole manifest, and cast Ua^v{q^ ''jP) in the form 
of Fig. ^ by setting 

(3.3) 



Ua^iv{q,r,p) =Ia{q) ] B^viq,r,p), 
where the nonperturbative quantity 



B^v{q,r,p) = Big^v + Biqtiqv + BiP^Pv + B^r^qv +Bsr^p^ 



(3.4) 
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+ 




Figure 3: The gray kernel (regular part with respect to q, and the composite massless excitation in the 
q'-channel. 



is the effective vertex describing the interaction between the massless excitation and two gluons. 
B^y{q,r,p) is to be identified with the "bound-state wave function" (or "BS wave function") of 
the two-gluon bound-state shown in Fig. ||, which, as we will see shortly, satisfies a homogeneous 
BSE. In addition, i/q^ is the propagator of the scalar massless excitation. Finally, Ia{q) is the 
(nonperturbative) transition amplitude introduced in Fig. 0, allowing the mixing between a gluon 
and the massless excitation; note that the imaginary factor "/" from the Feynman rule in Fig. ^ is 
absorbed into the definition of /ce(^). 
Evidently, by Lorentz invariance, 



and the scalar cofactor, to be referred to as the "transition function", is simply given by 



i{q) 



so that 



q- 



Vj{q,>;p) =I{q) ( ^ ]Bj{q,r,p)\ j= 1,...,5. 



(3.5) 



(3.6) 



(3.7) 



Note that, due to Bose symmetry with respect to the interchange jx ^ v and p ^ r,-we must 



have 



which implies that 



Biaiq,r,p) = -Bi^2iq,P,r).. 
Bi,2(0,-p,p)=0. 



(3.8) 
(3.9) 



4. Gluon mass and the BS wave-function: an exact relation 



The WI of Eq (2.11) furnishes an exact relation between the dynamical gluon mass, the tran- 
sition amplitude at zero momentum transfer, and the form factor B^. Specifically, contracting both 
sides of the WI with two transverse projectors, one obtains. 



P^'^{r)P^''{p)q"Va,v{q,r,p) = [m\r) - {p)]P^ {r)P'''\p) ■ 



(4.1) 



On the other hand, contracting the full expansion of the vertex ( |3.3[ ) by these transverse projectors 
and then contracting the result with the momentum of the background leg, we get 



q-p^^'^{r)P'''{p)Va^,{q,r,p)=il{q)[B,g^,+B2q^q,]P^'^{r)P'''{p), 



(4.2) 
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{A) 



Ia{q) 



V 



+ . 



■o- 




i ^ab 




Figure 4: (A) The vertex Uafiv is composed of three main ingredients: the transition ampHtude, la, which 
mixes the gluon with a massless excitation, the propagator of the massless excitation, and the (massless 
excitation)-(gluon)-(gluon) vertex. (B) The Feynman rules (with color factors included) for (i) the propaga- 
tor of the massless excitation and (ii) the "proper vertex function", or, "bound-state wave function", B^y 



where the relation of Eq (g/7|) has been used. Thus, equating both results, one arrives at 

iI{q)Bi{q,r,p) =m^{r) -m^{p) , B2{q,r,p) =0. 



(4.3) 



The above relations, together with those of Eq. ( |3.7| ), determine exactly the form factors Vi and V2 
of the vertex Va^v. namely 



V,{q,r,p)= ^ \ V2{q,r,p)=0. 

r 



(4.4) 



We will now carry out the Taylor expansion of both sides of Eq (4.3) in the limit ^ ^ 0. To 
that end, let consider the Taylor expansion of a function f{q,r,p) around ^ = (and r = —p). In 
general we have 

f{q, -p -q,p)= f{-p,p) + [2{q ■ p) + q']f'{-p,p) +2{q ■ pff"{-p,p) + ff{q^) , (4.5) 
where the prime denotes differentiation with respect to {p + qY and subsequently taking the limit 

Now, if the function is antisymmetric under p o r, as happens with the form factors Bi 2, then 
f{—p,p) = 0; thus, for the case of the form factors in question, the Taylor expansion is (/ =1,2) 

Bi{q, -p -q,p) = [2{q ■ p) + q'\B\{-p,p) + 2{q ■ pfB'l{-p,p) + ff{q^) . 



Using Eq (^^), and the corresponding expansion for the rhs, 

m^{r)-m^{p)=m^{q + p)-m^{p)=2{q-p)[m^{p)]' + ff{q^) 



(4.7) 



(4.8) 
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(a) 







(c) id) 
Figure 5: The complete BSE for the full three gluon vertex ir'„^^,(^,r,p). 



assuming that the 7(0) is finite, and equating the coefficients in front of {q- p), we arrive at 
(Minkowski space) 



[m\p)\ = imB\{p). 



(4.9) 



Note that this is an exact relation, whose derivation relies only on the WI and Bose-symmetry that 



V, 



,r,p) satisfies, as captured by Eq. ( |2.11[ ) and Eq. ( p^ , respectively. 



5. The Bethe-Salpeter equation 

As has become clear in the previous section, the existence of B[ is of paramount importance 
for the mass generation mechanism envisaged here; essentially, the question boils down to whether 
or not the dynamical formation of a massless bound-state excitation of the type postulated above is 
possible. As is well-known, in order to establish the existence of such a bound state one must (i) 
derive the appropriate BSE for the corresponding bound-state wave function, B^y, (or, in this case, 
its derivative), and (ii) find non- trivial solutions for this integral equation. 

The starting point is the BSE for the vertex Jr'f^^y{q,r,p), shown in Fig. ||. Note that, unlike 
the corresponding SDE of Fig. |2[ the vertices where the background gluon is entering (carrying 
momentum q) are now fully dressed. As a consequence, the corresponding multiparticle kernels 
appearing in Fig. ^are different from those of the SDE. 

The general methodology of how to isolate from the BSE shown in Fig. |5] the corresponding 
dynamical equation for the quantity B^y has been explained in [|l5[ |l8]]. Specifically, one separates 
on both sides of the BSE equation each vertex (black circle) into two parts, a "regular" part and 
another containing a pole l/q'^; this separation is shown schematically in Fig. ^ Then, omitting 
all other vertices, and the possible poles they too may have, the BSE for B^y{q,r,p) is obtained 
simply by equating the pole parts on both sides; specifically, [see Fig. ^ 



S-«= B"„'^AlP{k + q)AP-{k)J^, 



snmr 



(5.1) 



We will next approximate the four-gluon BS kernel by the lowest-order set of diagrams 
shown in Fig. ^, where the vertices are bare, while the internal gluon propagators are fully dressed. 
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{A) 






{B) 



77? , ^ 




k + a 



p + q 



771, jL 



Figure 6: (A) The separation of the vertex in regular and pole parts. (B) The BSE for the bound-state wave 
function B^y 



Going to Euclidean space, we define x = p^,y = k^, and z = {p + k)'^', then, after appropriate Taylor 
expansion, and use of the fact that 62 = [see Eq. (p3[)], the BSE becomes 



12712 



dyyB\{y)A'{y) 



dd sin^ 6 cos 6 



z+lO{x + y) + -{x^+y^ + lOxy) 
z 



Hz). 



(5.2) 



As a further simplification, we approximate the gluon propagator A(z) appearing in the BSE 
of (5.2) [but not the A2(};)] by its tree level value, that is, A(z) = 1 /z. Then, the angular integration 
may be carried out exactly, yielding 



25 jc 3 
4 y 4j2 



(5.3) 



6. Numerical analysis 



Next we discuss the numerical solutions for Eq. (5.3) for arbitrary values of x. Evidently, the 
main ingredient entering into its kernel is the nonperturbative gluon propagator, A(^). In order to 
explore the sensitivity of the solutions on the details of A(^), we will employ three infrared-finite 
forms, to be denoted by Ai(g), A2{q), and AT,{q), focusing on their differences in the intermediate 
and asymptotic regions of momenta. 

(i) Let us start with the simplest such propagator, namely a tree-level massive propagator of 
the form 

AY'{q^)=q^+ml, (6.1) 

where Mq is a hard mass, that will be treated as a free parameter. On the left panel of Fig. ||, the 
(blue) dotted curve represents Ai{q-^) for niQ = 376MeV. 
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p + q 





(01 ) 



Figure 7: The Feynman diagrams considered for the BS kernel. The interaction vertices are approximated 
by their tree level values, while the internal gluon propagators are fully dressed. 

(ii) The second model is an improved version of the first, where we introduce the renormalization- 
group logarithm next to the momentum q^, more specifically 



1 + - 



96712 



■In 



(6.2) 



where p is an adjustable parameter varying in the range of p € [2, 10]. Notice that the hard mass 
htq appearing in the argument of the perturbative logarithm acts as an infrared cutoff; so, instead 
of the logarithm diverging at the Landau pole, it saturates at a finite value. The (black) dashed line 
represents the Eq. (5.2) when p = 16, mo = 37 6 Me V, and = 4.3 GeV. 

(iii) The third model is simply a physically motivated fit for the gluon propagator determined 
by the large- volume lattice simulations of Ref. |j3|], and shown on the left panel of Fig. |[ The lattice 
data presented there correspond to a SU (3) quenched lattice simulation, where A(^) is renormalized 
at /X = 4.3 GeV. This gluon propagator can be accurately fitted by the expression 




(6.3) 



where mi{q^) is a running mass given by 



ml{q^) 



m 



q^ + p2rrfl 



(6.4) 



and the values of the fitting parameters are m = 520 MeV, gj = 5.68, pi = 8.55 and, p2 = 1.91. On 
the left panel of Fig. ||, the (red) continuous line represents the fit for the lattice gluon propagator 



■ O.U. 



given by Eq. ([6^31). Notice that, in all three cases, we have fixed the value of A^^ (0) 
Our main findings may be summarized as follows. 

fa) In Fig. ||, right panel, we show the solutions of Eq. (53) obtained using as input the three 



propagators shown on the left panel. For the simple massive propagator of Eq. (|6JJ), a solution for 
B[{q) is found for 01^ = 1.48; in the case of A2{q) given by Eq. (6^), a solution is obtained when 
a, = 0.667, while for the lattice propagator Ai,{q) of Eq. ( |6.3[ ) a non-trivial solution is found when 
a, = 0.492. 



(b) Note that, due to the fact that Eq. (5.3) is homogeneous and (effectively) linear, if B[{q) 
is a solution then the function cB[ (q) is also a solution, for any real constant c. Therefore, the 
solutions shown on the right panel of Fig. ^ corresponds to a representative case of a family of 
possible solutions, where the constant c was chosen such that B[ (0) = 1. 
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(c) Another interesting feature of the solutions of Eq. (5.3) is the dependence of the observed 
peak on the support of the gluon propagator in the intermediate region of momenta. Specifically, 
an increase of the support of the gluon propagator in the approximate range (0.3-1) GeV results in 
a more pronounced peak in B\ {q). 

(d) \n addition, observe that due to the presence of the perturbative logarithm in the expression 
for A2(^) and A3(^), the corresponding solutions B\ (q) fall off in the ultraviolet region much faster 
than those obtained using the simple Ai (q) of Eq. ( |6l] ). 

7. Conclusions 



In this presentation we have reported recent progress [ 27 ] on the study of the Schwinger mech- 
anism in QCD, which is the only self-consistent way to endow gluons with a dynamical mass. This 
mechanism relies on the existence of massless bound-state excitations, whose dynamical formation 
is controlled by a homogeneous BSE. As we have seen, under certain simplifying assumptions, 
this equation admits non-trivial solutions, thus furnishing additional support in favor of the specific 
mass generation mechanism described in a series of earher works [g, [T^, |ri[]. 

In the future it would be particularly important to consider the effects of bound-state poles 
in the SD kernels of not only the three-gluon vertex, as we did here, but of all other fundamental 
vertices of the theory. Such an investigation would eventually give rise to a coupled system of 
various homogeneous integral equations. Especially interesting in this context is the information 
that one might be able to obtain on the corresponding wave-function of the ghost-ghost channel. 
Specifically, according to the recent lattice findings [[T], Q ||, ^ ^], in the deep infrared the ghost 
dressing function F is finite, but the full ghost propagator diverges, a fact that strongly suggests 
that there is no dynamical mass associated with the ghost field (note that the finiteness of F can 
be easily accounted for by the presence of a gluon mass, saturating the perturbative logarithm of 



F [|10p). One would expect, therefore, that the solution of the corresponding system should give 



rise to a non- vanishing B\ , as before, but to a vanishing ghost-ghost wave function. 
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